Probability density functions of work and heat near the stochastic resonance of a 
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We study experimentally and theoretically the probability density functions of the injected and 
dissipated energy in a system of a colloidal particle trapped in a double well potential periodically 
modulated by an external perturbation. The work done by the external force and the dissipated 
energy are measured close to the stochastic resonance where the injected power is maximum. We 
Oh' show a good agreement between the probability density functions exactly computed from a Langevin 

dynamics and the measured ones. The probability density function of the work done on the particle 
satisfies the fluctuation theorem. 



I. INTRODUCTION 



^ ■ The study of fluctuations of the injected and dissipated power in a system driven out of equilibrium by an external 
1) [ force is nowadays a widely studied problem which is not yet completely understood. This is a very important 
^ ■ and general issue within the context of Fluctuation Theorems (FT) which constitute extremely useful relations for 
^ I ^ ' characterizing the probabilities of observing entropy production or consumption in out of equilibrium systems. These 
, relations were first observed in the simulations of a sheared fluid and later proved both for chaotic dynamical systems 
' and for stochastic dynamics]^. These works lead to different formulations which find powerful applications for 

^ ; , measuring free-energy difference in biology (see e.g. [p^ for a review) . The hypothesis and the extensions of fluctuation 
d • theorems jlTI have been tested in various experimental systems such as colloidal particles |l^, mechanical 
^ [ oscillators ]14|, electric circuits and optically driven single two- level systems [!£]. The effect of anharmonic potential 
I I on the motion of a colloidal particle has been tested by Blickle et al [|l3[ |l^, 17 1. In a recent experiment |^ it has 
' O ' been shown that FT holds for a colloidal particle confined in a double well potential and driven out of equilibrium 
near the stochastic resonance (see next section). The purpose of this article is to compare the probability density 
function (PDF) for work and heat measured in this experiment with those analytically computed from a non-linear 
Langevin dynamics. We also discuss the difference between the PDF obtained from a phase average of the driving and 
those obtained from a fixed phase only. Such a precise comparison between theory and experiment in a double well 
potential driven out of equilibrium has never been done before. There is only a numerical study, which has explored 
the distributions of the dissipated heat and of the work in a Langevin dynamics near the stochastic resonance |p^ , . 
fy-j I The paper is organized as follow. The properties of the stochastic resonance are recalled in the next section. The 
^\ . experimental set-up is described in section 3. The analytical PDF are derived in section 4. The comparison between 
■ the theoretical predictions and the experimental measurements is done in section 5. Finally we conclude in section 6. 

a^ 
o , 

QO n. STOCHASTIC RESONANCE 

o ■ 

A colloidal particle, confined in a double well potential, hops between the two wells at a rate r^, named the 
, , Kramers' rate, which is determined by the height SU of the energy barrier between the two wells, specifically = 
rS ] Tq^ exp(— ^^), where tq is a characteristic time, ks the Boltzmann constant and T the heat bath temperature [Q. 
. When the double well potential U is modulated by an external periodic perturbation whose frequency is close to 
rj. the system presents the stochastic resonance phenomenon^, i.e. the hops of the particle between the two wells 
synchronize with the external forcing. The stochastic resonance has been widely studied in many different systems 
and it has been shown to be a bona fide resonance looking at the resident time[^ ^, the Fourier transform of the 
signal for different noise intensity Q. Numerically, the stochastic resonance has been characterized by computing the 
injected work done by the external agent as a function of noise and frequency j^, j^]. In a recent experiment 
some of us have studied experimentally the Steady State Fluctuation Theorem (SSFT) in a system composed by a 
Brownian particle trapped in a double well potential periodically modulated by an external driving force. We have 
measured the energy injected into the system by the sinusoidal perturbation and we have analyzed the distributions 
of work and heat fluctuations. We find that although the dynamics of the system is strongly non-linear the SSFT 
holds for the work integrated on time intervals which are only a few periods of the driving force. In this paper we will 
compare this measured PDF with those derived analytically from a Langevin dynamics in which the exact potential 
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of the experiment has been used. 



III. EXPERIMENTAL SET-UP 



The experimental setup is composed by a custom buih inverted optical tweezers made of an oil-immersion objective 
(63x , N.A.=1.3) which focuses a laser beam (wavelength A — 1064 nm) to the diffraction limit for trapping glass beads 
(2 [m\ in diameter). The silica beads are dispersed in bidistillcd water in very small concentration. The suspension is 
introduced in the sample chamber of dimensions 0.25 x 10 x 10 mm'^, then a single bead is trapped and moved away 
from others. 

The position of the bead is tracked using a fast-camera with the resolution of 108 nm/pixcls which gives after treatment 
the position of the bead with an accuracy better than 20 nm. The trajectories of the bead are sampled at 50 Hz. 
The position of the trap can be easily displaced on the focal plane of the objective by deflecting the laser beam using 
an acousto-optic deflector (AOD). To construct the double well potential the laser is focused alternatively at two 
different positions at a rate of 5 kHz. The residence times (with i = 1,2) of the laser in each of the two positions 
determine the mean trapping strength felt by the trapped particle. Indeed if t\ = = 100/zs the typical diffusion 
length of the bead during this period is only 5 nm. As a consequence the bead feels an average double-well potential: 

Uq{;£) = ax^ — hj? — dx , (1) 

where a, b and d are determined by the laser intensity and by the distance of the two focal points. In our experiment the 
distance between the two spots is 1.45 /xm, which produces a trap whose minima arc at Xmin = ±0.45 ^m. The total 
intensity of the laser is 58 mW on the focal plane which corresponds to an inter-well barrier energy SUq = 1.8 fc^T. 
Starting from the static symmetric double-trap, (ti = T2) we modulate the depth of the wells at low frequency by 
modulating the residence times (r^) during which the spot remains in each position. We keep the total intensity of 
the laser constant in order to produce a more stable potential. The modulation of the average intensity is harmonic 
at frequency / and its amplitude {t2 — ti)/(t2 -I- ri), is 0.7 % of the average intensity in the static symmetric case. 
Thus the potential felt by the bead has the following profile in the axial direction: 

U{x,t) = Uoix) + Up{x,t) = Uq + c X sin(27r/t), (2) 

with axf^.^^ = 1.8 fc^T, hxl^^^ — 3.6 fcsT, d|xmi„| = 0.44 ksT and c|a;mi„| = 0.81 fc^T. The amplitude of the time 
dependent perturbation is synchronously acquired with the bead trajectory. The parameters given here are average 
parameters since the coefficients a, b and c , obtained from fitted steady distributions at given phases, vary with the 
phase {5a/ a w 10%, Sb/b « 5c/c w 5%). 

An example of the measured potential at t = jj and at < = ^ is shown on the Fig. |^a). The figure is obtained 
by measuring the two steady state probability distribution function P{x) of x, (corresponding to J7o + c and Uo — c 
respectively) and by taking U{x) = — In [-P(a;)]. 




FIG. 1: a) The perturbed potential ai t = jj and half a forcing period later, b) Example of trajectory of the glass bead and 
the corresponding perturbation at / = 0.1 Hz. 
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A. The equation of motion. 

The X position of the particle can be described by a Langevin equation: 

dU 

7i = -^+C, (3) 

with 7 = 1.61 10^^ N s the friction coefficient and ^ the stochastic force. The natural Kramers' rate (c 0) for 
the particle is = 0.3 Hz at T = 300 K. When c 7^ the particle can experience a stochastic resonance when the 
forcing frequency is close to the Kramers' rate. An example of the sinusoidal force with the corresponding position 
are shown on the figure Qd). 

B. The work and the heat. 

Since the synchronization is not perfect, sometimes the particle receives energy from the perturbation, sometimes 
the bead moves against the perturbation leading to a negative work on the system. 

In the following, all energies arc normalized by kgT. From the trajectories, we compute the stochastic Wg and the 
classical Wd works done by the perturbation on the system and the heat Q exchanged with the bath. These three 
quantities are defined by the following equations as in rcf . pc|] : 

Z"*"^*^ , dU(x,t) 

Wa{t,,tf) = - j'"^'' dtx ^^fj^^^ (4) 
f*°^'' , dU(x,t) 

where in this case i^- = j is a multiple of the forcing period. We use both Wg and Wd because they give complementary 
information on the fluctuations of the energy injected by the external perturbation into the system (see ref. and 
reference therein for a discussion on this point). For example, as discussed in refJ26l Ws/T is the total entropy 
production rate in this specific case ||2^. The heat and the work, defined in eq.|, are related through the first 
principle of thermodynamics: Q = ~AU + Wg, where AU = U{x{tf + to),^o + tf) ~ U{x{to),to); whereas the two 
works are related by a boundary term Wd = ~AUp + Ws, where AUp = Up{x{tf + to),tf + to) — Up{x{to),to). Since 
the characteristic time evolution of the perturbation is small compared to the fluctuation of position and due to the 
harmonic form of the perturbation, the integrals are computed as follows: 

tf/St 

Ws{to,tf) LO c 5t x{i) cos{ijj{to + tj)) 

i=l 

Wei{to,tf) = -AUp + W, (5) 
Qito,tf) = ~AU + Ws 

where St is the sampling time. We checked that the direct computing of integrals of Q and Wd gives the same results. 
It is important to stress that to can either take any value (as it has been done in rcf.|2g]) or be a multiple of 1//: the 
fluctuations in the two cases exhibit different PDFs, as we will see in the next sections. To compute the works and 
heat from experimental data for a given duration tf, we thus divide a single trajectory into different segments starting 
either with a fixed phase, or with different phases, before averaging the results over the whole trajectory, and then 
over different runs. In ref. |p6| the average work received over one period has been measured for different frequencies 
(tf = J in eq. ^). Each trajectory is here recorded during 3200 s in different consecutive runs, which corresponds to 
160 up to 7500 forcing periods, for the range of frequencies explored. It has been found that the maximum injected 
energy is around the frequency / « 0.1 Hz, which is comparable with half of the Kramers' rate of the fixed potential 
rK = 0.3 Hz. This maximum of transferred energy shows that the stochastic resonance for a Brownian particle is a 
bona fide resonance, as it was previously shown in experiments using resident time distributions ^ or directly 
in simulations ||]. It is worth noting that the average values of work in the case of a periodic forcing do not 
dependent on their definitions: only the boundary terms, which vanish in average with time, are different p^ . As the 
difference between Wd and Wg has been discussed in ref. we will focus here only on Wg that, in order to simplify 
the notation, will be indicated by W. 
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IV. EQUATIONS FOR THE W AND Q PDFS 

In this section we discuss the equation governing the time evolution of the work and heat PDFs. For a stochastic 
process described by eq. (H) the Fokkcr-Planck equation reads 

dtpix, t) = [U'ix, t)p] + ksTT^, (6) 

where p{x,t) is the PDF associated to the coordinate x, and F = I/7. Here and in the following, the prime denotes 
derivative with respect to x. 

Let us consider the joint probability distribution function of the position and of the stochastic work (j){x, W, t): in 
ref. it has been shown that the time evolution of such a function is governed by the partial differential equation 

with the starting condition 4>{x,W,to) = p{x,to)S{W). The unconstrained probability distribution of the work is 
given by 



= J dx<j){x,W,t). (8) 

By introducing the Fourier transform 

ipix,X,t) ^ J dWe-^^(j){x,W,t), (9) 

eq. becomes 

dtij{x, A, t) = T— [U'ix, t)^] + ksTT-^ - \—^, (10) 

with the starting condition ip{x,X,to) = p{x,to). Note that, by using the Fourier transform definition eq. (H), the 
wavenumber A associated to is a purely imaginary number, A = i|A|. 

Let us now consider the joint probability distribution ip(x,Q,t) of the position x and the heat Q exchanged by 
the brownian particle whose motion is described by eq. (||). The Fokker-Planck-like equation, governing the time 
evolution of such a function, reads p7| , 

dMx, Q, t) = (r(7V) - 9Q(rf/'V) - [{TU'kBT)dQv] - Sq [{VU'kBT)d,v] 

+rksT'^+TksTU-^,, (11) 

with the starting condition ip{x,Q,to) = p{x, to)<5(Q). Note that we use here the opposite sign convention for Q, eq. 
(^, with respect to that adopted in ref. |2^ . 

The unconstrained probability distribution functions of the heat reads 

<P{Q.t) = J dx^{x,Q,t). (12) 

Let x(a;. A, t) be the Fourier transform of ip{x, Q, t), as given by 

x(.T,A,i) = J dQe-^'^^ix.Q^t), (13) 

eq. ( pi] ) becomes 

dtxix, A, t) = ^ksT-^ + a, (Ft/'x) - XTU'^x - [(TU'kBT)x] - X(TU'kBT)d,x + A^FfesTt/'^, (14) 

with the starting condition x(x, A,<o) = Pix,to). In ref. p7| , it has been shown that by introducing the function 
g{x, A, t) defined by 



X(a;, A,i) = g{x,X,t)exp 



J±-^U{x,t) 



(15) 
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equation simplifies, as one gets rid of the first order derivatives with respect to x, and it becomes 

dtg = rksT^ - r/3^.g + ^U"g + ^—^g d,U, (16) 

with the starting condition as given by g{x, A, to) = p{x, to) exp[(/3 — 2X)U{x, to)/2]. 

We now show that a similar simplification can be obtained for eq. (|l^): let h{x, A, t) be defined as 

ipix, A, t) ^ h{x, A, t) cxp[~(3U{x, t)/2], (17) 

substituting h{x,X,t) into eq. (p^), we obtain the following equation for h{x,X,t): 

fj2 u TT/2 p — 9 A 

dth = rkBT^~Tp^h+-U"h+i^^hdtU. (18) 

The last equation is identical to (p^, however its starting condition reads h{x, A, to) = p{x, to) cxp[f3U {x, to)/2] , which 
is different from the starting condition of eq. (|l6|). 

As A is an imaginary number, one can split eq. ( |l6|) and eq. ([l8|), in a set of two equations for the real and imaginary 
part. For example, eq. (^8|) becomes 

dtliR = TkeT^ - rp^hn + ^U"hn + ( ^hn + \Mhi] dtU, (19) 



dthi = ThsT^ - rp^hi + ^U"hi + (l^hi - \X\hR ] dtU, (20) 



dx^ A ' 2 \2 

where /ir, and hi are the real and the imaginary part of /i, respectively. 

Equations (19) and (|20|), and the analogous equations for gn and gi, can be solved numerically for any value of 
|A|, and for any choice of the initial condition p{x, to) using, e.g^ MATHEMATICA The functions V'(a;, A, t) and 
x(a;, A,f), can be thus obtained by using equations ([l7| ) and (|l5|). The target functions 4>(x,W^t) and (p{x,Q,t) can 
then be obtained by taking the Fourier inverse transform, i.e. by inverting eqs. (^) and (p^, respectively. Also in this 
case the computation can be performed numerically. 



V. COMPARISON WITH EXPERIMENTS 



In this section, we compare the results for the work and heat PDFs, as obtained by solving the differential equations 
introduced in the previous section, with the experimental outcomes. We take the external driving frequency to be 
equal to / = 0.25 Hz, to have a good statistic, by allowing the observation of the system over a sufficient number 
of periods. Such a value is close to the natural Kramers' rate, which ensures that the system is in the stochastic 
resonance regime. 

We first consider the distribution of the work over a single period of the external potential (||) . The unconstrained 
probability distribution of the work ^{W,t), is obtained as follows. In order to solve eq. (1^, we solve numerically 
eq. ^) up to the time to ~ 50/ f ^ 1//, so as to ensure that the solution p(x,to) of eq. (^) represents the steady 
state distribution of the position of the particle. Such a solution is used as a starting condition for eq. (|lo|), which is 
solved numerically up to time ti = to + 1//, i.e. along a single period, and for different values of A, so as to obtain 
■0(a;, A,ti). The resulting unconstrained Fourier transform of the work PDF, defined as \1/(A, ti) = / da; i/;(a::, A, ti), is 
plotted in fig. |2| . As discussed in the previous section, by inverting eq. (^) and exploiting eq. (^, we finally obtain 
<^{W, ti). In fig. ^. the experimental histogram of the work done on the particle along a single period 1// is plotted; 
in the same figure the unconstrained probabifity distribution of the work ^[W,ti) is also plotted. We find a good 
agreement between the experimental distribution of the work and the expected PDF. It is worth noting first that 
we start sampling W only few periods 1// after the beginning of each experiment and second that the experimental 
results are averaged over different values of to in order to improve the statistics (either with a fixed phase or over 
different phases). 

We now check that the work PDF satisfies the SSFT, which we expect to be true for any time ti which is an integer 
multiple of the period 1// as discussed in jl^, ^ . Given the periodicity of the external driving, the stochastic work 
here considered corresponds to the total entropy production, as defined by Seifert in In figure^, we plot the 

symmetry function for the work PDF, defined as 



s[w, ti) = log [$(iy, ti)i^{-w, ti)] , 



(21) 
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FIG. 3: Diamonds: experimental probability distribution of the work done on the particle, along a single period, and with a 
single value of the initial phase. Full line: PDF of the work as obtained by numerical solution of eqs. ([7|)-(|lO[), with to = 50// 
and ti = to + 1/f. 



together with the corresponding work PDF, for two values of the final time ti. The fluctuation theorem states that 
the simmctry function has to satisfy the equality S{W, ti) = W, for any value of W. Inspection of the figure indicates 
that the fluctuation theorem is satisfied even at short times ti , the small deviations for the largest value of W being 
due to numerical and experimental difficulties in the sampling of the tails of the work PDF. More details on SSFT 
have been already discussed in ref. pof . 

We now consider the work distribution averaged over several phases. In order to obtain the probability distribution 
of the work, we solve eq. @ up to the time to{k) = (50 + k/N)/f, where = 20 and k = 0, . . .iV — 1, i.e. we 
consider A'^ initial conditions with a phase difference l/(iV/). For each value of k we solve eq. (|l0|), up to time 
ti{k) = to{k) + 1/f, so as to obtain ip{x, X,ti{k)\p{x,to{k)), i.e. the Fourier transform of the work PDF with the 
starting condition p{x,to{k)), and by antitransforming, we obtain 4>(x,W,ti{k)\p{x,to{k)). Finally the average work 
PDF is obtained as $(T4^) = l/N'^^kS dx(l}{x, X,ti{k)\p{x,to{k)). In fig. ||, this distribution is compared with the 
experimental outcomes. We find a good agreement between the expected PDF and the histogram of the work. 

Similarly to what has been done for the work, we now consider the PDF of the heat ^{Q), both with a single 
value of the initial phase, and averaged over several initial conditions. Also in this case, we use the function p{x, tg), 
solution of eq. (|^), as a starting condition for the numerical integration of eq. then by antitransforming the 

function x(a;,A,ti), we obtain ip{x, Q,ti) and finally ${Q,ti), as discussed in the previous section. The results 
for the case of a single value of the initial phase is plotted in fig. ||, where we show a nice agreement with the 
experimental data. In order to obtain the heat PDF averaged over several initial phases, we solve numerically eq. (|lj) 
in the time intervals [to(fc), ii(fc)], where to{k), and ti{k) have been defined above. Then by antitransforming the 
function x(a;, A,ti(/c)), we obtain (p{x,Q,ti\p{x,tQ{k)), and finally the unconstrained PDF of Q defined as ^(Q) = 
dx(f{x,Q,ti\p{x,to{k)). In fig. [7| we compare this expected PDF with the experimental outcomes, the 
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FIG. 4: Right panel: Symmetry function for the work PDF ^(W^, ti), eq. for two different values of t\. The symbols 

corresponds to the experimental data, the lines to the numerical solution of eqs. (|7[)-(|lo|) with = 50// . Red boxes and red 
dashed-dotted line: symmetry function obtained with = to + 2//. Blue diamonds and blue dashed line: symmetry function 
obtained with = to -f 4//. The full line corresponds to the expected behaviour S{W,t\) = W . Left panel: PDF of the work 
for ti = to -f 2// (red dashed-dotted line), and t\ — to + 4/ f (blue dashed line). 




FIG. 5: Diamonds: experimental probability distribution of the work done on the particle averaged over different initial phases. 
Full line, PDF of the work $(W) as obtained by numerical solution of eqs. (|7|)-([lo|), averaged over A'^ = 20 different phases, see 
text. 



agreement between the numerical PDF and the experimental data is good also in this case. 



VI. CONCLUSIONS 



In conclusion, we have experimentally and theoretically investigated the power injected in a bistable colloidal system 
by an external oscillating force. 

We have compared the PDFs of the stochastic work and of the heat measured in the experiment with those obtained 
theoretically. They exhibit large tails toward negative values and their shape are non-gaussian close to the resonance 
when averaging over a single period. 

In refs. [ p7| , it has been shown that if the external potential is quadratic, an analytical gaussian solution can be 
obtained for equation (|^) at any value of t and for cq. ( pi] ) in the long time limit. Here we have shown that the same 
equations can be solved numerically, also at finite time and for a more complex potentials, and they lead predictions 
on the system energy exchange that can be verified experimentally. Wc have shown that the Fourier transform of 
the heat and work PDFs obey the same differential equations (|l^) and (^8|), and the only difference lies in the initial 
conditions. This ensures $(W^) and ^(Q) to be different at short time, while one has <1>(VK) = <?(Q) in the long time 
regime, as one would expect, since (U) — for long times. 

Thus, our results suggest that the set of equations @) and (|l^) represents an useful and general tool to investigate 
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FIG. 6: Diamonds: experimental probability distribution of the heat exchanged by the particle, along a single period, and with 
a single value of the initial phase. Full line, ^(Q^ti) as obtained by numerical solution of eqs. (pl])-(p^), with to = 50// and 
ii = to + 1//. 
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FIG. 7: Diamonds: experimental probability distribution of the heat exchanged by the particle with the bath averaged over 
different initial phases. Full line, heat PDF ^'(Q) as obtained by numerical solution of eqs. (|ll[)-(|l4|), averaged over A'' = 20 
different phases, see text. 

the energy balance of microscopic systems in arbitrary potentials. 

Finally, we find that the work PDF satisfy the fluctuation theorem, even at short times. 
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